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Abstract 

The geometric measure, the logarithmic robustness and 
the relative entropy of entanglement are proved to be 
equal for a stabilizer quantum codeword. The entangle- 
■ i ment upper and lower bounds are determined with the 

■ generators of code. The entanglement of dual-containing 
| CSS codes, Gottesman codes and the related codes are 
. given. An iterative algorithm is developed to determine 

^ ' the exact value of the entanglement when the two bounds 
• are not equal. 

Index Terms: quantum code; Pauli measurement; 

■ multipartite entanglement 
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A quantum code encodes logical qubits in physical qubits. 
A quantum codeword is usually a multipartite state if 
the physical qubits are owned by different parties. The 
quantification of multipartite entanglement is basically 
open even for a pure multipartite state until now. How- 
ever, a variety of different entanglement measures have 
been proposed for multipartite setting. Among them are 
the (Global) Robustness of Entanglement pQ , the Rela- 
tive Entropy of Entanglement^ [3] , and the Geometric 
Measure [1] . The robustness measures the minimal noise 
(arbitrary state) that we need to add to make the state 
separable. The geometric measure is the distance of the 
state to its closest product state in terms of fidelity. The 
relative entropy of entanglement is a valid entanglement 
measure for a multipartite state, it is the relative entropy 
of the state to its closest separable state. The quantifica- 
tion of multipartite entanglement is usually very difficult 
as most measures are defined as the solutions to difficult 
variational problems. Even for pure multipartite states, 
the entanglement can only be obtained for some special 
scenarios. Fortunately, due to the inequality on the loga- 
rithmic robustness, relative entropy of entanglement and 
geometric measure of entanglement [5] [5] [7], these en- 
tanglement measures are all equal for stabilizer states 
[8]. A stabilizer state is a multiqubit pure state which 
is the unique simultaneous eigenvector of a complete set 
of commuting observables in the Pauli group, the latter 
consisting of all tensor products of Pauli matrices and 
the identity with an additional phase factor. 

A stabilizer state is the special case of a quantum code 
that encodes zero logical qubit. We may ask if the three 
entanglement measures are equal for a generic quantum 
codeword. The answer is true as shown in Section II. 
The three equal entanglement measures for a codeword 



then are simply called the entanglement of the codeword. 
The rest of the paper is organized as follows. In Section 
III, we derive the upper bound of the entanglement with 
Pauli measurements. In Section IV, the lower bound of 
the entanglement is obtained with the bipartition of the 
physical qubit system. Section V deals with the entangle- 
ment of the codewords of CSS codes. The entanglement 
of the family of Gottesman codes is the topic of Sec- 
tion VI. Section VII provides an iterative method for the 
possible exact value of the entanglement when the two 
bounds are not equal. Conclusions are drawn in Section 
VIII. 

2 The entanglement measure for 
quantum codewords 

The global robustness of entanglement R(p) PQ is defined 

as 

R(p)=mint (1) 
such that there exists a state A, satisfying 



a = (p + *A)/(l+t) € Sep 



(2) 



where Sep is the set of separable states. The logarithmic 
robustness is 



LR{p) =log 2 (l + R(p)). 



(3) 



The relative entropy of entanglement is defined as the 
"distance" to the closest separable state in terms of rel- 
ative entropy [3], 



E r (p) = min S (p\\ w) , 



(4) 



where S (p\\ ui) = —S(p)— tr{plog2co} is the relative en- 
tropy, S(p) is the von Neumann entropy.The geometric 
measure of entanglement for pure state is defined as 



Eg(\i>)) = 



1 0)6 Pro 



io g2 i(0i v>r 



(5) 



where Pro is the set of product states. An exten- 
sion of the definition for mixed state p is also available, 
E g (p) = min^gsep — log 2 tr(pui), however, E g is an en- 
tanglement monotone only for pure states p = \?p} . 
It has been shown that the maximal number N of pure 
states in the set {\ipi) \i = 1, -/V}, that can be discrim- 
inated perfectly by LOCC is bounded by the amount of 
entanglement they contain[B]: 



log 2 N<n- LR(\iPi)) <n- £ r (|Vi)) <n~ E g {\^)), 



(G) 



where n — log 2 Dh , Dh is the total dimension of the 
Hilbert space, and x = jj Yli=i x i denotes the "average" . 

An n-qubit stabilizer state \S) is defined as a simul- 
taneous eigenvector with eigenvalue 1 of n commuting 
and independent Pauli group elements Mj. The n eigen- 
value equations M, \S) = \S) define the state \S) com- 
pletely (up to an arbitrary phase). The group generated 
by the product of the n operators Mj is called the sta- 
bilizer S, and Mi are generators of S. A subgroup of S 
with n — k generators is also called a stabilizer [12] . de- 
noted as M. However, M C S stabilizes a 2 k dimensional 
space. In principle, such a space is the coding space 
{\tp) , s.t.T \ip) = , VT e M}, corresponds to a sta- 
bilizer code encoding k into n qubits. In addition to the 
n — k stabilizer generators, a stabilizer code also has log- 
ical operations X\, . . . , Xk and Z\, . . . , Zk- We can take 
the basis codewords for this code to be 

|0) = l[T\0f n , 

TeM 

|c> = xf.-.X^Io), 

where c =(ci, . . . , C&) is a binary vector. Zi |0) = |0) for 
i = 1, . . . , k. 

For a stabilizer state \S), it has been shown that [8 

LR(\S)) = E r (\S)) = E g (\S)). (7) 

This is also true for a stabilizer quantum codeword |c) . 

Proposition 1 The entanglement measures of the loga- 
rithmic robustness, the relative entropy of entanglement 
and the geometric measure of entanglement are all equal 
for a stabilizer quantum codeword |c) . 

Proof: It is enough to prove the statement for |o) , 
since |c) is locally equivalent to |0) . Notice that |o) is 
stabilized by {Mi, . . . , M n _fe, Z\,..., Zf.}- The n — k gen- 
erators and k logical Z operators commute with each 
other and are independent. Thus |0) is a stabilizer state, 
the three entanglement measures are equal according to 
Ref. H]. 

For a stabilizer quantum codeword, we will simply call 
these three quantities the entanglement, and denote as 
E(\c)). 

3 Entanglement upper bound 

A generator which is composed of Z operator or identity 
in each qubit and does not contain X or Y operators in 
any qubit is called a Z-type generator. 

Proposition 2 The entanglement of a codeword is up- 
per bounded by the minimal number of stabilizer genera- 
tors which are not Z-type generators. 

Proof: The codeword |0) = N]J^(I + 
Mi) \0)® n , where N is the normalization factor. For a 
Z-type generator M h we have (7 + M;) |0)® n = 2 |0)®" . 
In the operator product Yl^I + Mi), we may move 



the factor (I + Mi) to the rightmost. The number of 
terms R in the linear decompositions of |0) into product 
states is upper bounded by 2 r , with r' the number of 
non Z-type generators. Let r = minr' be the minimal 
number of non Z-type generators. Notice that we can 
artificially increase the number of non Z-type generators 
by replacing a Z-type generator with the product of the 
Z-type generator and a non Z-type generator. Hence we 
count the minimal number of non Z-type generators. It 
follows that the Schmidt measure [S] 

E s — min log 2 R 

is upper bounded by r. The geometric measure is upper 
bounded by Schmidt measure [TU], hence the theorem 
follows. 

For a stabilizer group M, up to over all phase 
±l,±i, each generator Mi = X ai Z bi , with X ai = 
(gtj I°' J , Z bi = (g) . Z^ i] , where Oi and &j are the bi- 
nary vectors (a fl , a i2) . . . , a in ) and (bu, b i2 , ...,b in ), re- 
spectively. An alternative representation of the generator 
Mi is (a^l bi) . We may use generator matrix ( follow Ref. 
PTj . it is called stabilizer matrix in Ref. [T^] ) (A\ B) 
to represent the stabilizer group M, where A and B are 
(n — k) x n matrices with elements Ay = ay, Bij = bij.lt 
is always possible to arrange (A\B) in the form of (see 
e.g. [H] Ch.4) 





D 


F 




(o 





J 





(8) 



by the permutations of the qubits and replacements of the 
generators with other elements in the stabilizer group. 
Here I is an r x r identity matrix, with r the F 2 rank 
of A, where F 2 denotes the integer field {0, 1} with ad- 
dition and multiplication modulo 2. With the stan- 
dard form ([SJ of the generator matrix, we may improve 
the upper bound of the entanglement. We then inves- 
tigate the effect of Pauli measurements on codewords. 
We work on the generators of M in the first r lines 
of (j8|) and neglect the other n — k — r Z-type gener- 
ators. We have M x = X ® M[ or M t = Y ® M{, 
M l+1 =Z® M' t+l oi M <+ i = I ® M' i+1 (i = 1, . „, r - 1) 
for the standard form of the generators. Denote |0„_i) = 

N'TliZiV + Mi+i) |0)® (n_1) ,where N' is the normaliza- 
tion factor, then Pauli Z measurement on the first qubit 
will project the codeword |0) to 

P { z +\0) - |0>®|0 n _i), (9) 
P Z W |0) = |l)®M{|0 n _i>. (10) 

The projection operators on the j th qubit are P^2 = 
1(7 ± Zj). The two measurement results ±1 are equally 
probable. Similarly, the X or Y measurements on the 
first qubit also project the codeword to two equally prob- 
able states corresponding to the two measurement re- 
sults ±1, except for possible special case of X or Y mea- 
surements with only one result. When M\ = X ® M[, 
the Pauli X, Y measurements on the first qubit will 
project the codeword |0) to |0) = ±(|0) ± |1)) <X>(/± 



M[) |O n _i) , Pjg |0) = §(|0) ±i\l)) ®{lTiM[) |O n _i) . 
When Mi = Y (8 M( , the Pauli X, Y measurements on 
the first qubit will project the codeword |0) to P^± |0) 
= i(|0>±|l>] ®{I±iM[) |0 n _x) , P£2 |0) = j(|0>±i |1» 
®(/ ± M[) |0 n _i) . There is the case that |0„_i) is the 
eigenvector of M[, so occasionally one of the two out- 
comes is annihilated. 

As a property of the Schmidt measure [§] , for any se- 
quence of local projective measurements that finally com- 
pletely disentangles the state vector in each of the 
measurement results, we obtain the upper bound 



ea ) i 



(ii) 



where N mea is the number of measurement results with 
non-zero probability. 

The minimal number of local Pauli measurements to 
disentangle a stabilizer quantum codeword can be called 
its " Pauli persistency" . 

Proposition 3 The entanglement is upper bounded by " 
Pauli persistency" for a stabilizer quantum codeword |c) . 

Proof: With the formula (TTT1) . the fact that different 
measurement results of codeword are obtained with prob- 
ability 1/2, and geometric measure is upper bounded by 
Schmidt measure, the statement follows. 

Example. — Consider [[8,1,3]] code [T3] with stabilizer 
in the standard form 
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the Pauli Z measures are applied at the 1,5,7 qubits. 
For each qubit measured, the corresponding row and col- 
umn are deleted. What left for the remain qubits is the 
stabilizer 

X I z z z 
I X Z I z 
Z Z X I i • 
Z I I X I 



(12) 



The | O5 ) generated by stabilizer (TT21 is a product of graph 
state |G 4 ) with |0) , |0 5 ) = |G 4 ) ® |0) . The graph state 
\Gi) is generated by a new stabilizer obtained with delet- 
ing the last column of (fl2)) . "Pauli persistency" of graph 
state IG4) is 2, so "Pauli persistency" for codeword of 
[[8, 1, 3]] code in Grassl code-table (T3] is 5. 

We will use proposition [5] to obtain the entanglement 
upper bound of codeword in the following except Table 
1, where proposition [3] is used for tighter upper bounds. 

4 Entanglement lower bound 

The index of physical qubits is denoted as I = 
{1,2, ... ,n}, for a bipartition, we may assign m qubits 



to A, and the remain n — m qubits to B. The index 
sets of A and B are I4 and = I— Z4, respectively. 
The reduced state of the codeword |0) then should be 
Pb = Trj( |0) (0| . The bipartite entanglement of the 
bipartition{Z.4,2g} then will be — Trpg log 2 pb, the en- 
tropy of p B . 

Proposition 4 The entanglement of a codeword is lower 
bounded by any bipartite entanglement of the codeword, 



E > -Trp B \og 2 p B - 



(13) 



Proof: If we define E r t,i as the relative entropy of en- 
tanglement with respect to some bipartition, we have 
that E r > E r bi since the set of fully separable states 
is a subset of the bipartite separable states. Notice that 
E r = E and E r i,i is equal to the bipartite entanglement 
Ebi = —TrpB log 2 ps for pure state, the statement then 
follows. 

We will obtain the entropy of ps by diagonalizing pb 
and at last the entropy can be expressed with the code 
stabilizer. The entanglement of the codeword is lower 
bounded by the maximal bipartite entanglement among 
all bipartitions. 

Since a Z-type generator does not contribute new items 
to codeword |0) , we simply ignore Z-type generators. 
Thus we take {A \B) = (I D\E F) in the following. A, B 
are rxn binary matrices with r the number of non Z-type 
generators and r < n — k. 

The codeword 



iV^(-l) tt ^X^|0)®" 



X 



jeiA 



iei B 



where the summation on r dimensional binary vector jU is 
from (0, 0, ... , 0) to (1, 1, ... 1), (pA)j is the j th compo- 
nent of the binary vector fiA, N is the normalization fac- 
tor, and a(p) — 'Yl li< i{p-iai){P'ibJ) . We my rewrite a(p) 
as 

a( M ) = i[/^> T - Tr(ATA)} = ^pT^ , (14) 

with A = diag{pi, . . . , p r }, and Ti is the matrix T with 
diagonal elements nullified, where 



T = AB = F 1 



DG 



The convention for binary addition is mod 2. T is 
symmetric for any two generators should commute with 
each other, namely, AB T + B T A = 0. The reduced 
state pb = E^n^Wc/W-l)^ ^ 



ing normalization, which is 



i(n — m) y(fJ,'A) 



X}, disregard- 



pb = EII U- l (^) ) -(- 1 ) aW+a( "' ) IHe)(H B |, 



with \(jiA)g) = UpA) 2 



(Mk) 



We then consider to diagonalize in order to ob- 
tain its entropy. Without loss of generality, let I\ = 1, 
T 2 = 2, I m = m < r, and denote /i = (v, r),with 
v = (//!,..., /i m ),r = (/i m+ i,...,ju r ). then|(^4) e ) = 
(M-D)i, • • • , G"-D)n-r) = |t, ^£>) • Denote 
l*M> = E T (-l) a(Al) l(M)B),then 

= E ^,'(-l) Q(Al)+Q(M,) |(M)6>((M) B | 

ZAl/' ,T,t' 



For the orthogonality of ^(z/)) , we turn to 

= Z)(-l) O( " ,T)+0,(,/ ' T) ^,r)iJ,(^,r)iJ. (15) 
r 

The bipartite entanglement of the codeword is at least m 
when all | v I / (^)) are orthogonal with each other. This is 
obviously from the factor that 



1 



(1,1,. .,1) 

y 

i/=(0,0,...,0) 



where |4 r (t')) are orthonormal and the normalization fac- 
tor is retrieved. The bipartite entanglement may be less 
than m only when the non-orthogonality of the ensemble 
^(f)) is found. The conditions for nonzero (^(v 1 ) l^i^)) 
are 



(v + v',0)D = 0, 
(y + v')T 5 = 0, 



(16) 
(17) 



where stands for the r — m dimensional zero vector 
(0,0,..., 0). T 3 is produced by deleting the first m 
columns and the last r — m rows of the r x r matrix Ti , 
so r 3 is a m x (r — m) submatrix of Ti. More explicitly, 
we may write Ti as 



1 3 



r 3 
r 4 



Then a(i/,r) + a(i/,r) = \(yY 2 v T + v'Y 2 v' T ) + (v + 
v')Tj,t t + tT^t t . The T4 term always contributes a +1 
factor in the summation of Eq. (ll5[) for T4 is symmetric 
and with nullified diagonal elements. The T 2 term con- 
tributes a constant factor in the summation of Ea. ((T5|) . 
Then Eq. flTTl) follows. Let D' be the matrix produced 
by deleting the last r — m rows and preserving the first 
m rows of D, the rank of the m x (n — m) matrix 
Q(A, B) — (r 3 , D') gives the number of independent vec- 
tors Y&(v)) for a specific bipartition of first m qubits for 
A with respect to last n — m qubits for B. Maximizing 
with respect to all bipartitions except the m > r cases, 
we hence obtain the maximal of bipartite entanglement 
as the lower bound of the entanglement 



Since Q(A, 6)isamx(n- m) matrix, its rank must not 
exceed min{m, n — m} < I 2J 7 so we have E\ < I §1 . 

One of the special case that should be notified is when 
m = r. It follows that r 3 is an r x matrix and does not 
exist at all. So Eq. pTl) disappears, and we only need to 
consider Ea. lflTj)) . If rank^^D — r, then Ea. lfTo]) fulfills 
only when v = v' , hence Ei = r. 

5 Entanglement of CSS codes 

5.1 Dual-containing CSS codes 

An important class of quantum codes, constructed from 
classical codes, invented by Calderbank, Shor [T5] & 
Steane [IS], has the generator matrix of the form (e.g. 



(A\B) = 



U 





V 



(19) 



Ei = max rankw 2 Q(A,B), 

•partitions 



(18) 



where U and V are I x n matrices. Requiring UV T = 
ensures that the generators commute with each other. 
As there are 21 stabilizer conditions applying to n qubit 
states, k = n — 2l qubits are encoded in n qubits. We may 
write the classical parity check matrix U in a systematical 
way 

U = [ I D ]. (20) 

Since AB T = 0, we get (x(fi) = for all binary vectors 
fi. Consider the case of m = I, we have v = fi. Then 
(^(f') |*(^)> = 5 V D,u'D- From vD = u'D, we have (v + 
v')D = 0. Thus the condition for the orthogonality of 

vD = => i> = 0, (21) 

for all binary vector v. The lower bound of the entan- 
glement of the codeword is / when the condition ([2T]) is 
fulfilled. For dual-containing code, we have V = U, thus 
UU T = 0, so that DD T = I, the condition flUJ) is ful- 
filled. The lower bound of entanglement is Ei = I = 

The upper bound of the entanglement E u of the code- 
word |0) is the number of X generators now, which is I. 
Thus the entanglement of dual-containing CSS codeword 
is 

n — k , „ 

■ (22) 

For a CSS code that is not dual-containing, the upper 
bound E u is still Z, the number of X generators. The 
lower bound is the binary rank of D. 

5.2 The graph state of a CSS code 

A stabilizer code with stabilizer generators 
Mi, . . . , M„_fc and logical operations Xi,...,Xk 
and Z±, . . . , Zk, is equivalent to the codeword stabi- 
lizer (CWS) code [T7] defined by codeword stabilizer 
{Mi, . . . , M n -k, Zi, . . . , Zk} and word operators which 
are products of Xj. Any CWS code is locally Clifford- 
equivalent to a standard form of CWS code with a 
graph-state stabilizer and word operators consisting 
only of Z operators. The standard codeword stabilizer 
is generated by XiZ Vi . The set of YiS forms the 



adjacency matrix of the graph [17) . Hence, given a 
quantum stabilizer error-correcting code, we can always 
find the corresponding graph state. The entanglement 
of the codeword of the quantum stabilizer code and the 
graph state should be the same, since they are locally 
Clifford-equivalent. A CSS code has a generator matrix 
(fTTH) and U can further written in the form of (|2"U)) . We 
now construct the graph state stabilizer. The generator 
matrix of {Mi, . . . , M„_ fe , Z t ,..., Z k } is 



(23) 



where (0 \W) is the generator matrix for logical op- 
erations Zi,...,Zk. With elementary row transforma- 
tion we have transformed U into the systematical form 
\ I D 1 . We want show that it is always possible to 

transform Y Tr into the form of \ D' I 1 , where I 
W J L J 

is an (n — I) x (n — I) identity matrix. We first trans- 

r V 1 . r i 
form into [it r J , where P is an upper trian- 




gle square matrix, namely Py 



case that Pjj — 0, we then interchange the j 



for i > j. There is the 
th qubit 

with some later qubit such that Pjj = 1. This is always 
possible since the elements of the j — th line of P can 
not be all zeros, otherwise the elements of j — th line of 
[ R P ] should be all zeros due to mutual commuta- 
tion of the generators, namely, UV T = and UW T = 0. 
That is 



IR 1 



DP T = 0. 



(24) 

From which we can deduce that if the j — th line of P are 
all zeros, we have Rji = 0, for all i < I. An all zero line in 
[ R P 1 means the generator is the identity, this is not 
the case. It is easy to transform [ R P ] to [ D' I ] , 
then Eq. (O will be 

D' = D T . 

Performing Hadamard transformation to the last n — I 
qubits, the generator matrix undergoes the transforma- 
tion 





D 






















D T 


!)■ 




I 


D T 


o) 



Thus the adjacency matrix of the locally Clifford- 
equivalent graph state of CSS codeword is 



7 



D 
D T 



(25) 



A graph state with adjacency matrix of (|2"5j) is two- 
colorable, we can simply assign the first I qubits with 
one color and the remain n — I qubits with another. The 
entanglement upper bound should be [TS E u = n— (n — 

fcf 



I) = I. The lower bipartite bound is [5] E\ 
rankp 2 D. 

For a dual-containing CSS code, the corresponding 
graph state is further characterized by DD T = I in addi- 
tion to two-colorable. The D matrix has a full rank and 



the upper and lower bounds of entanglement coincide. 
The entanglement of the codeword is also E = I with the 
theory of graph state and formula (|25l) . 

5.3 Toric Codes 

Toric code is proposed to encode quantum information in 
topological structure [19], it is a kind of quantum LDPC 
code [TB] . The toric code is based on a k x k square lattice 
on the torus. Each edge of the lattice is attached with a 
qubit, so there are n = 2fc 2 qubits. For each vertex s and 
each face p, operators of the following form are defined: 



A s — ]Q Xj, Bp 

j£star(s) 



n 



j ^boundary (p) 



These operators commute with each other. Due to fj A s 
= 1 and Yip Bp = 1, there are m = 2k 2 — 2 independent 
operators constitute the stabilizer of the toric code. The 
code encodes n — m = 2 qubits. Toric code is a kind 
of CSS code from its definition. We will show that it is 
not a dual containing code, the upper and lower bound 
of entanglement may not coincide. By proper numbering 
the edges, the generator matrix can be written in the 
following form 



U = 



I 

I I 



i n 



n 



i i 
r r 



v = 



n 1 



I 

I I 



n 1 



Tl 



I I 

V V 



where I is the k x k identity matrix, Q is a k x k matrix 
with 2k nonzero entries, 



i 

i i 



i 



i 



i 



where I 1 , fi' and J7 T 'are the matrices obtained by delet- 
ing the last line of I, f2 and f2 T , respectively. In order to 
show that the D matrix has not a full rank in general, we 
transform it into the following form by elementary row 
transformation: 



A 



A A 1 



where 1 =(1, 1, 
vector, A = [/&_ 



. . , 1) T is a k — 1 dimensional column 
L , 1]. The (fc 2 - 1) x (k 2 + 1) matrix D 



is apparently not a full rank matrix, for the rank of f2 is 
k — 1, we have the lower bound of entanglement 



1. 



While the upper bound of entanglement is E u — fc 2 — 1. 
Hence the upper bound is no longer equal to the lower 
bound unless k = 2. 

The codeword of a toric code is a highly entangled state 
for large k, the entanglement scales as E ~ ~. It seems 
that the area law [3U] does not work for toric code. Area 
law usually means that the entanglement is proportional 
to the boundary area when the bulk of qubits is cut to 
two parts. The number of bulk qubits now is n = 2k 2 , 
so according to area law, the (bipartite) entanglement 
should be proportional to the cutting length L, which 
is now proportional to k in most cases. However, the 
largest bipartite entanglement is proportional to fc 2 , cor- 
responding to a sophisticated cutting curve with length 
L cx k 2 , despite that a random cutting of toric qubits 
into two parts usually yields a boundary length L oc ^/n,. 
Hence the area law still works, but the boundary may be 
very long. So care should be taken when we talk about 
the area law of multipartite entanglement and the largest 
bipartite entanglement. 

6 Entanglement of Gottesman 
codes and the related codes 

6.1 Gottesman codes 

A serial quantum codes [[2 m , 2 m -m-2, 3]] (m > 3) that 
fulfill quantum Hamming bound had been proposed by 
Gottesman [21] . By construction, the first two generators 
of the stabilizer are X\ ■ ■ ■ X 2 ™ and Zy ■ ■ Z 2 ™ . An ex- 
plicit construction of the remaining m generators is given 
by the matrix (H\CH), where H = [h , hi, /i2 m -i] 
with the (fc + l) th column hk being the binary vector 
representing integer fc (fc = 0,1,..., 2 m — 1) and C is 
any invertible and fixed point free m x m matrix, i.e., 
Cs ^ and Cs ^ s for all s e F™ except s = 0. The 
generator Z\ ■ ■ ■ Z 2 ™ is a Z-type generator and omitted 
hereafter regarding the entanglement of codewords. We 
may arrange H — [Hq,Hi, H 2 , ■ . ■ H m ] , with Hq = ho = 



[0,0,. 



,or, #1 = 

is an m x 



\ho. m -i,ho 



,h 2 ,hi] 



In 



matrix whose column vector has 



and // , 

3 r 

weight j. The generator matrix of the last m generators 
in the standard form will be ( / D | F G ) , with 
D=[ H 2 H 3 ■■■ H m ],F = C,G = CD, 



F 1 



DG 



0. 



(26) 



where we have used the facts that any two rows of matrix 
H are orthogonal, and each row vector of H has even 
weight, so J2T=i HiHT = 0, and £™ 2 H t Hf = H X H^ = 
1. Eq. (|2"6"|) and the fact that the generator Ai...A2>™ 
does not contain any Z operator leads to 

a{jj) = 0, 



for all m + 1 dimensional binary vectors /1. To obtain 
the lower bound of entanglement for Gottesman code- 
words, we should verify if condition (|2"Tj) is satisfied or 
not. When the generator Xi ■ ■ ■ X 2 m is considered, the 
whole D matrix for m + 1 generators is 



D 



1 

H 2 H3 



(m - 1) F2 



where 1 and are vectors (1,1,..., 1) and (0, 0, ... , 0) 
with proper dimensions, respectively. We have 



DD 1 



Em' (m\ Y^m' -. „T 

i=l \ 2 i) l^i=l LI1 2i 



with 



[m/2\ . Notice that the I element of the m 



dimensional vector H 2 il T is the weight of I th line of H 2 i, 
each line of H 2 i has the same weight ti by the definition 
of H, each column of H 2 i has the same weight 2i, so 
mti — 2z(™) is the total weight of the matrix H 2 i, thus 

U = (Szi) , and H 2l l T = ET=! U = ET=i = 

2 m - 2 , which is in F 2 for m > 3, so that YT=i = 

(F. Meanwhile ^ti Q = Y,£o GD " 1 = ^ " ^ 
which is 1 in F2 for m > 2. We have 

DD T = I. 

The condition (|2"T]) is fulfilled. The entanglement lower 
bound of Gottesman codeword is m + 1. The number of 
the generators which are not Z-type is m+1, so the upper 
bound of the codeword is m + 1. We conclude that the 
entanglement of codewords is m + 1 for Gottesman code 
[[2™, 2 m - m - 2, 3]] (m > 3). Written with the length of 
the code n = 2 m , the entanglement of the codewords is 



E = log 2 n + 1. 



6.2 Family of 8m codes 



(27) 



The family of codes with parameters [[8m, 8m — l m — 5, 3]] 
with l m = |~log 2 m] was constructed [2"2"]. Alternative gen- 
erator matrices of the codes were given [33] based on 
Gottesman codes. The number of generators is l m + 5. 
There is one Z-type generator in the stabilizer. Thus 
the upper bound of entanglement (might not be tight) 
of codewords is E u = l m + 4. To obtain the lower bound 
of entanglement, we will utilize the generator matrices of 
[2"5] directly instead of transforming them into the stan- 
dard form of Eq.©. The code can be divided into m 
blocks, each block has 8 qubits. The generator matrices 
can be written as 

(A\B) = (A 1 ,A 2 ,...,A m \B 1 ,B 2 ,...,B m ) (28) 

where Ai and Bi are (l m + 5) x 8 binary matrices, and 
every line of {Ai \Bi) is either a line from the gener- 
ator matrix of Gottesman [[8,3,3]] code or corresponds 
to 



It is observed that AjAj 



0,AiBj = Q,B,Bj = 0, for all i,j . Thus we have 

T = AB T = A i B J = 0- ( 29 ) 



Hence = for arbitrary binary vector fi. Meanwhile, 
we have J2i A%Af = 0, thus 



A A = 0. 



(30) 



Notice that elementary row transformation of A keeps 
Eg. (1301) . After elementary row transformation, A can be 
transformed to standard form 



Ah A' 



I D 




(31) 



Hence 



and 



A' A' 



I + DD T 




0. 



DD 1 



I. 



(32) 



What is crucial is the dimension of the identity matrix in 
Eq.([3T]). There is one obvious Z-type generator. Besides 
this one, it is always possible to work out an full rank 
identity matrix from A by elementary row transforma- 
tion, thus the rank of the identity matrix in Eq. (pnT) as 
well as in Eq.fl32]) is l m +4. Based on Eq.® and Eci. (|52|) . 
the lower bound of the entanglement is Ei = l m + ^ which 
coincides with the upper bound. Hence the entanglement 
of the codewords of length n = 8m code is 



E= \\og 2 n\ +1. 



(33) 



Notice that Eq. (|2"T)l for the entanglement of Gottesman 
codes can be merged into Eq. (|3"3"|) . 

6.3 Pasted codes 

The [[13, 7, 3]] code is obtained [23] by pasting Gottes- 
man [[8,3,3]] code and cyclic [[5,1,3]] code. For the en- 
tanglement of the codewords of [[13, 7, 3]] code, there are 
one Z-type generator among the 6 generators, the upper 
bound should be E u — 5. A direct calculation shows that 
the lower bound Ei is also 5. So we have the entanglement 
E = 5, which fulfill Eq.([3"3l- 

The family of perfect codes [[n m ,n m — 2m, 3]] with 
n m — (4 m — 1)/3 for to > 3 is successively constructed by 
pasting Gottesman 2 2 ( m_1 ' code (occasionally, we denote 
the code with its length when it is not confusing) with 
7i m _i code [23] [25] . For the entanglement of codewords 
of n m code, there are 2m generators and one of them is 
Z-type generator, the upper bound is not difficult found 
to be E u = 2m — 1. For the entanglement lower bound, 
let's first consider [[21,15,3]] code obtained by pasting 
Gottesman 2 4 code with cyclic [[5,1,3]] code. We may 
utilize the special case mentioned at the end of Section [3J 
then we have Ei = rankf 2 D = 5, notice that the Z-type 
generator has already removed. Similarly, for n m code, we 
have Ei = 2m — 1. The entanglement is E — 2m — 1 and 
can be written with respect to the code length n = n m 
as 



E = [log 2 n] 



(34) 



Another family of codes [[8n m , 8n m — 2m — 3,3]] with 
to > 2 is successively constructed by pasting Gottesman 



2 2m+1 code with 8n m _i code [25] • The entanglement 
upper bound should be E u ~ 2m + 2 which is determined 
by the number of non Z-type generators of Gottesman 
2 2m+1 code according to the structure of 8n m _i code. 
The entanglement lower bound is Ei = 2m + 2. Since 
the code is a pasting of several Gottesman codes, the 
entanglement lower bound can be obtained to be Ei — 
2m + 2 by a technic similar to that of Subsection l6.2l The 
entanglement is E = 2m + 2 and can also be written in 
the form of Eq. (|3"4f with the code length n = 8n m . 



7 Iteration Algorithm 

Denote / = (0| &s), the closest product state \&s) = 
®j{ x ] 1°) +Vj I 1 )) witn \zj\ 2 + \y 3 \ 2 = 1- Using Lagrange 
multiplier method, we have L = |/| 2 — Ylj ^j(\ x j\ + 
\yj\ 2 — 1), where Xj are the multipliers. The extremal 
equations should be §£-J* - X 3 x* = 0, §^f* - \ jV * = 0. 
Let Zj — Uj/xj, we have 



df/dyj 
df/dxj 



(35) 



Thus the group element M^Mg 2 



■ M^_ k k is isomor- 



phic to (j2i=i J2i=i /"A) = (pA\ (J.B) , where 

= (/ii, /^2, ■ • ■ , /in-fc) is the binary vector. \iA and fj,B 
are binary vectors of length n. So that 



i — k 



N< 



U(I + M t )®(x 3 \0)+ yj \l)) 



3 = 1 



f^(o| to z" fl i" A (-i) a WH)'" s 

/x=0 



)(^|0)+ % |1» 



= A^(-l) tt ^H)^(0|® n A^ A 

/i=0 



)(Xj 



= A^(-l) a ^(-z)^ ff Yl 



„i-(M)i„(M)3 

J i Hi 



.(36) 



where g — (gi, . . . ,p n _fe), and gi is the number of Y 
operator in Mi. From (|35[) . the iteration equation for Zj 



EH(M) J= i(-i) a(Ai) ("0"- g n m#J -^ 



(37) 



Notice that the iteration may sometimes fail to reach the 
global maximum of |/| 2 . So, if the ultimate iteration 
result of the separable state |<I>s) is the closest product 
state for 1 0) , the entanglement of a quantum code from 



Table 1 The entanglement and the bounds 



[[n,k,d]] 


E 


E u 




[[4,1,2]] 


2 


2 


2 


[[4,2,2]] 


2 


2 


2 


[[5,1,3]] 


2.9275 


3 


2 


[[5,2,2]] 


2 


2 


2 


[[6,1,3]] 


2.9275 


3 


2 


[[6,2,2]] 


3 


3 


3 


[[6,3,2]] 


2 


2 


2 


[[6,4,2]] 


2 


2 


2 


[[7,1,3]] 


3 


3 


3 


[[7,2,2]] 


4 


4 


3 


[[7,3,2]] 


4 


4 


3 


[[7,4,2]] 


3 


3 


3 


[[8,1,3]] 


5 


5 


4 


[[8,2,3]] 


4.8549 


5 


4 


[[8,3,3]] 


5 


5 


4 


[[8,4,2]] 


4 


4 


4 


[[8,5,2]] 


3 


3 


3 


[[8,6,2]] 


2 


2 


2 


[[9,1,3]] 


5 


5 


4 


[[9,2,3]] 


5 


5 


4 


[[9,3,3]] 


5 


5 


4 


[[9,4,2]] 


4 


4 


4 


[[9,5,2]] 


3 


3 


3 


[[9,6,2]] 


2 


2 


2 



iteration method will be 

2 



E = - log 2 |/„ | = n - k - n s 
l 

-21og 2 '" 



E (-irc>n*J 



where rt s is the number of Z-type generators, /*, a;j» and 
are the extremal values of /, Xj and yj, respectively. 
The entanglement for some quantum codes with listed 
generators by Grassl 13 is as Table 1. The calculation is 
based on the iterative algorithm for unequal upper and 
lower bounds of entanglement. 



8 Conclusions 

The three entanglement measures (the geometric mea- 
sure, the logarithmic robustness and the relative entropy 
of entanglement) are proved to be equal for quantum sta- 
bilizer codeword. The entanglement upper bound of a 
stabilizer codeword can be the minimal number of non Z- 
type generators. A Z-type generator is the tensor product 
of identity and/or Pauli Z operators. Further tight upper 
bound is the " Pauli persistency" , the minimal number of 
Pauli measurements to resolve the entanglement. The 
entanglement lower bound based on bipartite entangle- 
ment is reduced to a formula of calculating the , maximal 
rank of some matrices. The matrices are derived from 
the generator matrix of the code stabilizer. The entan- 
glement of a self-dual CSS code is proved to be the num- 
ber of X generators regardless of the detail structure of 



the stabilizer. We also derive the adjacency matrix of 
the corresponding graph state of CSS code. Upper and 
lower bounds of entanglement are given for toric codes, 
the entanglement is about half of the code length. Com- 
ments are given on the area law of entanglement for toric 
codes. The entanglement values of the Gottesman codes 
[[2 m ,2 m — m — 2,3]] (m > 3), 8m codes and Gottesman 
pasting codes are equal to their minimal numbers of non 
Z-type generators. The entanglement E of the Gottes- 
man codes and the related codes scales with the code 
length n as E = [log 2 n\ + 1 or E — [log 2 n\ + 1. An 
iterative algorithm is developed to obtain the entangle- 
ment of the codeword as precisely as possible. 
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